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ABSTRACT 

Hydrodynamical winds from a spherical two-temperature plasma sur- 
rounding a compact object are constructed. The mass-loss rate is computed 
as a function of electron temperature, optical depth and luminosity of the 
sphere, the values of which can be constrained by the fitting of the spectral 
energy distributions for known X-ray binary systems. The sensitive depen- 
dence of the mass loss rate with these parameters leads to the identification 
of two distinct regions in the parameter space separating wind-dominated 
from non wind dominated systems. A critical optical depth (tc), as a func- 
tion of luminosity and electron temperature, is defined which differentiates 
these two regions. Systems with optical depths significantly smaller than Tc 
are wind-dominated. 

The results are applied to black hole candidate X-ray binary systems in 
the hard spectral state (Cyg X-1, GX 339-4 and Nova Muscae), and it is 
found that the inferred optical depth (r) is similar to suggesting that they 
are wind regulated systems. On the other hand, for X-ray binary systems 
containing a neutron star (e.g., Cyg X-2) r is much larger than Tc indicating 
the absence of significant hydrodynamical winds. 

Key words: binaries: accretion, accretion-discs - black hole physics - hydro- 
dynamics 

1 INTRODUCTION 

X-ray binary systems containing a black hole candidate are typically found to be in two 
different spectral states. In the hard state, the broad band X-ray spectrum can be described 
as a power-law (photon spectra index ^ 1.5) with a high energy cutoff around 100 keV. 
In the soft state, the spectrum consists of two components. There is usually an extended 
power-law (with spectral index ^ 2.5) and a soft X-ray emission which has a spectral shape 
similar to that of a black body. For a recent review of the observations and phenomenological 
description of these sources see Tanaka & Shibazaki (1996). 

The modeling of the hard state spectra can be described in terms of an unsaturated 
Comptonization process of soft photons in a region with hot electrons {kT ?a 50 keV) and 
electron scattering optical depths of order unity. In a pioneering study, Shapiro, Lightman 

6 Eardley (1976) identified this hot region with a geometrically thick, optically thin, hot 
accretion disc. In this model, the gravitational energy dissipated in the disc, heats the 
ions which in turn transfer the energy to electrons by Coulomb interactions. However, the 
electron-ion Coulomb interaction rate is inefficient in such an environment of low density 
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and high electron temperature, and this leads to a large difference in temperature between 
the electrons and ions, with the ion temperature reaching nearly virial values {VS^^K). The 
importance of radial advection of energy in such a disc was noted by Ichimaru (1977). Taking 
advection into account Narayan & Yi (1994) constructed self-similar solutions for the disc 
equations called Advection Dominated Accreting Flows (ADAF) showing not only that the 
proton temperatures approach their virial values, but also that the radiative efficiency of 
accretion can be significantly reduced as a result of the advection of energy into the black 
hole. In an alternative description Chakrabarti & Titarchuk (1995) argued that under certain 
conditions a shock may arise in such accretion discs and identified the hot Comptonizing 
region with the post-shock flow. Despite the differences in the geometry, radiative processes 
and the detailed disc structure, both these models have in common the presence of a two- 
temperature plasma. 

Such a plasma is a natural outcome of any accretion disc model which (a) identifies the 
hard state X-ray spectrum as a result of the unsaturated Comptonization process of soft 
photons, (b) assumes that the viscous energy dissipated heats the ions and (c) that the only 
mechanism for energy transfer between the ions and electrons is Coulomb interaction. We 
note here that these assumptions may not be valid since the viscous energy dissipated may 
heat the electrons preferentially if a strong equipartition magnetic field is present in the disc 
(Bisnovatyi-Kogan & Lovelace 2000). In addition, there could be unknown mechanisms which 
transfer energy between ions and electrons more efficiently than the Coulomb interaction. 
Thus, it will be useful to have an independent observational signature which could confirm 
the existence of two- temperature plasmas in black hole candidate systems. 

The nearly virial proton temperature of this plasma suggests the possibility of a strong 
hydrodynamical wind arising from these systems (e.g., Piran 1977; Takahara, Rosner & 
Kusunose 1989; Kusonose 1991). Such an outflow could transport away a signiflcant fraction 
of mass, energy and/or angular momentum, thereby affecting the structure and stability of 
the disc and affecting the radiative efficiency of accretion for a given mass transfer rate. 
Chakrabarti (1999) and Das (1999) have studied the possibihty of outfiows in the context of 
the shock/centrifugal barrier models. They find that for certain values of the disc parameters 
(e.g., accretion rate and specific entropy) a hydrodynamical wind occurs. For the ADAF 
disc solutions, Blandford & Begelman (1999) argued that only a small fraction (< 1%) 
of the gas actually falls into the black hole and the rest is driven away as a wind (in 
an advection dominated inflow-outflow solution - ADIOS), thereby effectively reducing the 
radiative efficiency. Bcckcrt (2000) confirmed this result for different viscosity laws while 
Quataert & Narayan (1999) showed that the X-ray spectra from such a wind driven accretion 
process can explain the observed spectra of some black hole candidate systems in quiescence. 
These calculations were undertaken for low mass accretion rates, and it is not clear how the 
system will behave in the high accretion rate regime inferred for black hole candidate systems 
in the hard state. 

The formation of hydrodynamical winds from accretion discs depends on the structure 
and geometry of the discs. Thus, detailed calculations of the outflow are intrinsically model 
dependent. Further, rehable calculations of the structure of such discs are difficult because 
uncertainties exist in the vertical distribution of energy dissipation in the disc associated 
with our lack of detailed understanding of the viscosity. We note that internal magnetic fields 
in the disc could also facilitate (or inhibit) the formation of winds. Here, electromagnetic 
forces may accelerate and coUimate the wind to form high velocity jets as observed in mi- 
croquasar systems (Mirabel & Rodriguez 1999). Considering these uncertainties a prudent 
approach would be to estimate the wind characteristics using only those parameters which 
can be directly constrained by the fitting of spectral energy distributions. Such an analysis 
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would allow a rough estimation of the magnitude of the mass and energy lost in the form of 
a wind for a system in question. With such an objective in mind, we report in this paper on 
calculations of the structure of an hydrodynamical wind similar to those found by Takahara 
et al (1989), but in the context of a uniform spherical two-temperature plasma around an 
accreting black hole. Rather than treating the entire disc/wind configuration with detailed 
heating and cooling processes as in Kusunose (1991) the calculations are parameterized in 
terms of the electron temperature of the cloud, Tg, the optical depth, r, and the total lumi- 
nosity of the source. We adopt this approach as these parameters are constrained by spectral 
fitting analyses. In the next section, the formulation of the problem and the numerical re- 
sults are presented. The application of these results to observed systems is given in §3 and 
discussed in the last section. 



2 HYDRODYNAMICAL WINDS 

The mass loss rate due to a hydrodynamical wind depends sensitively on the ion temperature. 
This allows a system to be defined as wind-dominated if the ion temperature is greater than 
some critical value (T^), while for temperature less than this value the mass and energy loss 
rate is small enough not to affect the dynamics of the system. This critical temperature 
will scale as the virial one i.e. kTf, oc GM/R where R is the size of the region and M is 
the mass of the black hole. For a two-temperature plasma the power output is mediated by 
electron-proton Coulomb interactions and the luminosity, L, is given by 

L/V = ^n{kTi - kTe)Uep (1) 

where V = A/SnR'^ is the volume of the sphere, n is the average number density, Tj and 
are the ion and electron temperatures, and z/gp is the frequency of the electron-ion Coulomb 
interaction. A critical optical depth Tc{= narR) may be defined for the system, wherein for 
systems with r < Tc, the ion temperature is larger than Tc (see below). Taking the above 
relationships into account with the Coulomb exchange rate given as (Spitzer 1962) 

= 2.4 X lO^HnApT-^/'^ (2) 

where lnA{^ 15) is the Coulomb logarithm and p is the mass density, it follows that 

OC LV2r3/4^-l/2 (3) 

Since L and are parameters which can be constrained by spectral fitting, Tc can be directly 
estimated for a black hole of a given mass. Spectral fitting can also give information about 
the optical depth of the system, which can then be compared with Tc, to determine whether 
the system could be wind-dominated or not. The above analysis shows that Tc does not 
depend on the size of the system R. This is fortunate since R is not well constrained by 
observations. 

To quantify the above analysis we solve for the hydrodynamical wind structure from a 
spherical two-temperature plasma with an optical depth (r), electron temperature Tg and 
luminosity L. The central mass is taken to M = lOM© and the radius of the sphere is fixed 
ai R — 20GM/c^. We use the basic equations of the hydrodynamical theory of stellar winds 
(Parker 1958) which are the conservation of radial momentum, 

dv _ _dP _ CM 

^ dr d,r ^ 

and the conservation of mass, i.e. the mass outfiow 

Mo — pvA'Kr'^ — constant. (5) 
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Here v is the radial velocity, P = Kp^ is the pressure and F is the adiabatic index. These 
equations can be combined to give, 

dr- |„-£l] 

where Cg = {dP / dpY^"^ is the sound speed. At the sonic point, both the numerator and de- 
nominator vanish. Wind solutions were constructed by integrating from the sonic point. The 
inner boundary is taken at R where the number density n — rj [Rot) and ion temperature is 
calculated using equation (1). The solution is constructed such that at infinity both the den- 
sity and pressure tend to zero. A wind solution docs not exist if the ion temperature at the 
surface k% < kTmin ~ {^)GMmp/R. For high values of kT^ > kT^ax ~ {^)GMmp/R, 
the flow is supersonic at the surface i.e. the sonic point is actually located at a radius less than 
R. Thus we restrict this analysis only to those values of F such that kTmin < kTi < kTmax- 
Note that for F = 5/3, kTmin ~ kTmax and there does not exist a hydrodynamical wind solu- 
tion with a sonic radius greater than R for any kT^. Thus this analysis is restricted to values 
of F < 5/3. Such regimes have also been studied by Chakrabarti (1999) and for example, 
in the limit of efficient conduction of heat an isothermal wind flow results corresponding to 
F ?^ 1. Radiative processes can also lead to redistrubution of energy in the wind leading to 
F < 5/3 (Chakrabarti 1999). 

In Figure 1, the calculated mass loss rate as a function of ion temperature is illustrated 
for a flxed optical depth (r — 1). The virial temperature T^ = GMrUp/kR for the sphere 
is 5 X lO^-*^ K. As expected, Mg is sensitive to the ion temperature: a factor of two 
increase in Tj causes an increase of M„ by two orders of magnitude. This justifies defining 
a critical temperature beyond which the system is wind-dominated and below which it is 
not. For comparison, we note that the ion temperatures in ADAF type solutions correspond 
to Tj ~ T^, which in the presence of winds are lowered (Misra & Taam 2000). The mass 
loss rate can be expressed in terms of luminosity, optical depth, and electron temperature. 
For a typical luminosity (L = 10"^^ ergs/s) and electron temperature (Tg = 50 kcV), the 
variation of mass loss rate with optical depth is shown in Figure 2. Here, the sensitivity of 
Mo with T is apparent: a factor of two increase in r decreases Mg by at least three orders of 
magnitude. We show the contour plots of Mg for the r, Tg and r, L planes in Figures 3 and 
4 respectively. These plots highlight the steep variation of with all three parameters. 

The r, Tg and L parameter space can then be divided into two regions corresponding 
to high and low values of M^. In particular we define the wind dominated region as corre- 
sponding to Mo > Merit s-nd the complementary region as the non wind dominated one. The 
choice of Merit is ad-hoc and is chosen to be 10^^ g/s. Since the dependence of Mg on the 
parameters is steep, choosing Merit to be lO^'' g/s still leads to a division of the parameters 
space to nearly similar regions. A critical optical depth Tc can be defined as a function of 
L and T^ for which Mg — 10^^ g/s. The topology of Figures 3 and 4 and the qualitative 
analysis described in the beginning of the section (eqn 3) suggests that Tc can be represented 
as: 

r = A( ^^'^ ^ (7) 

^50kcV^ 4038crgs/s^ ^ ' 

We constrain A, a and (3 for different values of Merit, T, M, and R and present the result 
in Table 1. 

From these results it can be seen that A, a and /5 do not vary greatly with changes in 
Merit, r, M or R. Thus to within a factor of two, the critical optical depth can be represented 
as 
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Table 1. Values of A, f3 and a in equation (|7|). Here is in units of g/s and R is in units of = IGMjc? . 



log Merit r 


M(M0) 


Rirs) 


A 




/3 


18 


1.05 


10 


10 


1.9 


0.90 


0.60 


17 


1.05 


10 


10 


2.5 


0.84 


0.56 


18 


1.2 


10 


10 


1.7 


0.85 


0.56 


18 


1.05 


30 


10 


1.3 


0.86 


0.57 


18 


1.05 


10 


100 


3.7 


0.90 


0.58 


rc = 2( 




jO.85/ 


L 


•jO.55 






•SOkeV 


' 4038 


ergs/s 









(8) 



We can now use equation (H) to determine from observations if a particular system is wind 
dominated or not. 

3 APPLICATION TO X-RAY BINARIES 

Assuming that a two-temperature plasma configuration is established subject to the con- 
ditions described in §2, we discuss the possible applications of hydrodynamical winds to 
binary systems containing compact objects. In the following we first discuss the results for 
the black hole candidate systems Cyg X-1, GX339-4, and Nova Muscae. For convenience, 
the results are summarized in Table 2. 

The hard state spectrum of Cyg X-1 has been fitted by Gierlinski et al. (1997) with 
a Comptonization model, and they obtained the following parameters: r ^ l,kTf. ^ 140 
keV and L = 3.5 x 10^'' ergs/s (for a distance of 2.5 kpc). Using these values we find from 
equation (^ for a black hole of IOMq that Tc = 2.6 which is on the order of the observed 
r 1. Similarly for another black hole system, GX 339-4 in the low state, Zdziarski (1998) 
obtained r = 1.93, kTe = 48 keV and L = 3 x 10^^ ergs/s. For these values Tc = 3.6 which is 
again of the order of the observed r. Thus, our analysis of both Cygnus X-1 and GX 339-4 
in their hard states indicates that these systems are on the borderline of having a strong or 
weak wind. This may point to the existence of a feedback mechanism wherein the energy 
and mass loss from the wind regulate the hot disc structure. 

Unlike the persistent systems discussed above, black hole X-ray novae are transients. In 
these sources, the luminosity increases by several orders of magnitude in 10 days and 
then decays exponentially in a time-scale of ~ 1 month. During the peak of the outburst 
the system is usually in the soft spectral state and makes a transition to the hard state 
in ^ 2 — 3 months (see Tanaka & Shibazaki 1996). Observations of the X-ray transient 
source Nova Muscae (GS 1124-68) by the Ginga satellite through most of its evolution were 
reported in Ebisawa et al. (1994). Since the Ginga energy band is restricted to < 20 keV, 
the roll-over of the spectrum at high energies was not observed. Thus, a Comptonization 
fit does not constrain both and r. In particular only a low energy spectral index (a^) 
(during different times of the evolution) was measured. For a Comptonized spectrum, is 
approximately related to Te and r by, 

3 /9 4 

-(3 + a,) = --±^- + - (9) 

where, y = {4kTe/meC^) max(r, r^) is the Compton y-parameter. For Nova Muscae in the 
hard state, was observed to ~ 0.5 and at a luminosity, L = 4.3 x 10'^^ ergs/s (Ebisawa et 
al. 1994). If we now assume that T^, ^ 100 keV one obtains from equation (|^), r ~ 1.7. The 
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Table 2. Comparison of inferred and critical optical depth for various X-ray binaries. The luminosity, electron temperature 
(Te) and optical depth (r) have been quoted from the references. The symbols WD,WR and NW stand for wind dominated, 
wind regulated and no wind respectively. For GS 1124-68 (Nova Muscae) and Cygnus X-1 (soft state) the electron temperature 
is assumed and the optical depth has been estimated using equation (^. Note that for the soft state, the assumptions used in 
this analysis may not be valid (see text). 



Source Name 


Luminosity (ergs/s) 


fcTe (keV) 


r 




Comment 


ref. 


Black Holes 


Hard State 


Cygnus X-1 


3.5 X lO^'' 


140 


1 


2.6 


WR 


1 


GX 339-4 


3 X 10^'' 


48 


1.9 


3.6 


WR 


2 


GS 1124-68 


4.3 X 10^^ 


100 


1.7 


0.6 


NW 


3 


GS 1124-68 


4.3 X 10^'^ 


200 


1.2 


1.1 


WR 


3 


Soft State 


Cygnus X-1 


2.9 X lO^'' 


200 


0.40 


3.3 


WD 


4 


GS 1124-68 


2.5 X 10^^ 


200 


0.37 


3.0 


WD 


3 


Neutron Stars 


Low luminosity 


XB 1608-52 


4.1 X 10^'^ 


7 


7.6 


0.13 


NW 


5 


XB 1636-536 


2.7 X 10^^ 


1.8 


16.7 


0.12 


NW 


5 


XB 0748-676 


3.7 X lO^e 


2.4 


24.1 


0.04 


NW 


5 


XB 1254-69 


7.8 X 10="^ 


1.9 


15.6 


0.06 


NW 


5 


XB 1820-30 


1.4 X 10^^ 


3.5 


13.2 


0.15 


NW 


5 


High luminosity 


XB 1820-30 


4.7 X lO^* 


3.3 


11.9 


1.2 


NW 


5 


Cyg X-2 


1 X 10^* 


3.7 


9.4 


0.5 


NW 


5 


GX 17-1-2 


2.5 X 10^8 


3 


13 


0.8 


NW 


5 


GX 9+1 


1.7 X 10^8 


3 


11 


0.6 


NW 


5 


GX 349-1-2 


2.2 X 10^8 


3.7 


10 


0.9 


NW 


5 



References; 1: Gierlinski et al. (1997), 2: Zdziarski et al. (1998), 3: Ebisawa et al. (1994), 4: Gierlinski et al. (1999), 5: White,Stella 
& Parmar (1988). 

critical optical depth for such a system is then Tc = 0.6 which is less than the inferred r, 
indicating the absence of a strong wind. However, a strong wind is indicated if one assumes 
kTe ~ 200 keV instead, since the inferred r ^ 1.2 which is close to the critical value of 
Tc = 1.1. Thus, hydrodynamical winds can be important in the hard state of Nova Muscae 
provided that the electron temperatures are ^ 150 keV. 

In the soft state of black hole candidate systems the power-law spectrum is steeper than 
the hard state. Gierlinski et al. (1997) analyzed the soft state Cygnus X-1 data and found 
that the energy spectral index ~ 1.5 and the power-law extends up to ~ 200 keV with no 
apparent cutoff. They argued that thermal Comptonization does not describe the spectra 
well. Instead, the power-law is probably due to non-thermal Comptonization. In this case, the 
presence of non-thermal electrons indicates that the electrons are probably heated directly 
instead of mediated by processes involving protons. Thus, the basic assumption underlying 
this study is probably not valid for the soft spectral state of Cygnus X-1. Nevertheless, from 
the observed = 1.4, L = 3 x 10^'' ergs/s and assuming that kTe ~ 200 keV, one would 
infer using equation @ that r 0.4. This is much smaller than the critical value of Tc = 3.3, 
perhaps indicating the presence of a strong wind. A similar result is obtained from the soft 
state of Nova Muscae (Table 2). 

In addition to black hole candidates, the spectra from X-ray binary systems containing 
a neutron star can also be described as being due to Comptonization. Hot two temperature 
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accretion flows around neutron stars have been recently constructed by Medvedev & Narayan 
(2000). On the other hand, the hot region may not be the accretion disc itself but could be 
instead an extended corona surrounding the boundary layer between the disc and the surface 
of the star. This analysis can be still applied to these systems as long as the assumption 
that the total luminosity of the source is channeled by the protons to the electrons remains 
valid. The critical optical depth when M = IAMq and R= 10 kms is found to be 

50 keV lO"** ergs/s 

Spectral fitting of EXOSAT data from several neutron star binaries has been undertaken 
by White, Stella & Parmar (1988) which are summarized in Table 2. For Cygnus X-2, 
they constrain kTe ~ 3.7 keV, r = 9.4 and L = 10^^ ergs/s (for a distance of 8 kpc). Using 
equation ( [I0| ) we find that Tc = 0.5 which is significantly less than the observed value. Similar 
results were obtained for other neutron star systems in both in the high and low luminosity 
levels (Table 2). This implies that unlike black hole candidate X-ray binary systems, their 
neutron star counterparts are not dominated by a hydrodynamic wind. This result is further 
supported by the analysis undertaken by Medvedev & Narayan (2000) who constructed 
self-similar accretion flows onto neutron stars and found that hydrodynamic winds are not 
important in these systems. 



4 SUMMARY AND DISCUSSION 

The possible occurrence of hydrodynamical winds in X-ray binary systems has been investi- 
gated. From simple considerations, the conditions under which these winds can be important 
have been identified. It is found that the mass loss rate in such winds depends sensitively 
on the luminosity of the source as well as the electron temperature and optical depth of the 
coronal region. The steep dependence of the mass loss rate on these parameters facilitates 
the use of a critical optical depth to indicate whether a given system can support such a 
wind. Application of the theory indicates that winds can exist in the hard state of the black 
hole candidates Cyg X-1, GX 339-4 and GS 1124-68. 

Strong winds in these systems may decrease the radiative efficiency {rj = L/MiC^, where 
Mi is the mass inflow rate) by carrying away a substantial amount of matter and energy. For 
Nova Muscae, the radiative efficiency has been estimated to he r] = 0.01 for the soft state 
and rj = 0.05 for the hard state (Misra 1999). These values are lower than that expected 
from an Keplerian disc ( ~ 0.1) which could be due to the presence of a strong wind rather 
than energy advection into the black hole. Winds may also effect the thermal stability of 
accretion discs by introducing an additional channel for energy loss (e.g., Piran 1977). 

We reiterate that the analysis undertaken in this paper is based on the assumption that 
the gravitational energy dissipated in the system is transferred to the electrons by the ions 
via Coulomb interactions. This naturally restricts the analysis to only certain systems where 
this is valid. As mentioned earlier, the soft state spectrum of black hole candidate systems 
is probably of a non-thermal origin and, hence, the analysis performed in this study cannot 
be applied. Even for the hard spectral states of compact X-ray binaries considered here, 
alternate models to a two -temperature plasma description have been proposed. For example, 
the X-ray spectra could be due to magnetic flare activity above a cold disc (Poutanen & 
Fabian 1999) or produced by a disc with a rapidly varying radial temperature profile (Misra, 
Chitnis & Melia 1998). Thus, it should be emphasized that the results presented here are 
specifically discussed within the framework of a two-temperature plasma model. 

The spherical geometry assumed here is simplistic even though two-temperature discs 
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are geometrically thick. Further the effect of angular momentum or of convection on the 
dynamics of accretion and the wind outflow have not been taken into account. The former 
effect can increase the mass outflow rate as a result of centrifugal support, whereas the 
latter effect can affect the thermal structure of the underlying disk and, hence, the existence 
of a wind. In this context, the thermal structure is dependent on the direction in which 
angular momentum is radially transported by convection and the magnitude of the viscosity 
parameter (see Narayan, Igumenshchev, & Abramowicz 2000). A more detailed analysis 
should take these effects and radiative heating/cooling of the wind into account in at least 
two spatial dimensions (e.g., Chakrabarti & Molteni 1993). However, such analyses will, by 
necessity, be model dependent and limited by the uncertainties in the vertical structure and 
geometry of the hot disc. 

The hydrodynamical winds described in this paper may be confined to form a jct-like 
structure by the geometry of the disc and magnetic fields. Further if they are accelerated to 
relativistic speeds by electromagnetic forces, they may provide the origin of the radio jets 
observed in the black hole candidate systems known as microquasars (Mirabel & Rodriguez 
1998). On the other hand, the radio jets may be a different phenomenon unrelated to hy- 
drodynamical winds. In that case it is desirable to have direct observational signatures of 
these outflows. Since the Thomson optical depth of the winds calculated here (i.e. the optical 
depth from the surface of the sphere R to infinity) is typically < 0.1, column densities of the 
order of 10^^ cm~^ are indicated. Although the X-ray continuum spectra are not expected 
to be altered by the outfiow, a significant fraction of the column density may not be highly 
ionised giving rise to observable absorption and/or emission lines. The detection of P-Cygni 
type profiles by high resolution X-ray satellites Chandra and Newton-XMM may provide 
evidence for the existence of such winds and serve as a useful diagnostic not only of the 
wind structure, but also the physics underlying the disc as well. 

RM acknowledges support from the Lindheimer Fellowship at Northwestern University. 
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Figure 1. The variation of the mass outflow rate with ion temperature for adiabatic indices V = 1.01 (nearly isothermal) and 
r = 1.2. The other parameters are: M = WMq, R = 20GM/c'^ and t = tiutR = 1. 
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Figure 2. The variation of the mass outflow rate with optical depth for adiabatic indices F = 1.01 (nearly isothermal) and 
r = 1.2. The other parameters are: M = IOMq, R = 20GM/c^, kTe = 50 keV and L = 10^* ergs/s 
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Figure 3 
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Figure 3. Contours of mass outflow rate (Mo) for luminosity and optical depth. The contours correspond to from top to 
bottom Mo = 3x IQi*, 10i**,3 X lO^'^, lO^'^, 3 x IQi^, IQi^ g/s. The other parameters are: M = IOMq, R = 20GM/c^, F = 1.01, 
kT, = 50 keV. 
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Figure 4 
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Figure 4. Same as in Figure 3, except that the contours are for electron temperature and optical depth for a luminosity of 
L = 103* ergs/s. 
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